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Abstract 

Nobeyama Milhmeter Array was used to observe milhmeter-wave afterglow of 
GRB 030329 at 93 GHz and 141 GHz from 2003 April 6 (8.23 days after the burst) to 
2003 May 30 (61.97 days). A sensitive search for CO (J = 1 — 0) emission/absorption 
from the host galaxy of GRB 030329 was also carried out. Unresolved millimeter 
continuum emission at the position of GRB 030329 was detected until 2003 April 21. 
We found a steep decline of continuum flux (oct^^ °) during this period, in accord with 
a previous report. Moreover, our data implies that the decay was accompanied by 
possible plateaus phases, or bumps, on a time scale of several days. From an integrated 
spectrum, produced by summing up the data from 2003 April 10 to 2003 May 30, we 
found a possible emission feature, which could be a redshifted CO {J = 1 — 0) line. 
Its position and redshift coincide well with those of GRB 030329, though further 
observations are required to confirm the detection. If the emission feature is real, the 
observed CO flux is 1.4±0.52 Jy km s~^, corresponding to a large molecular gas mass 
of M(H2) > 10^ Mq. This implies that the host galaxy, which is optically faint, is 
highly obscured due to a rich interstellar medium. 
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1. Introduction 

It is now widely believed that long-duration gamma-ray bursts (GRBs) are extremely 
energetic phenomena associated with star formation in galaxies at a cosmological distance. 
Detection of the supernova component SN 2003dh (Stanek et al. 2003; Hjorth et al. 2003; 
Kawabata et al. 2003) in the afterglow of the HETE-2 localized GRB 030329 (Vanderspek et 
al. 2004) has established the association of some of the GRBs with core-collapse supernovae 
from massive stars. GRBs are often hosted by dusty, star-forming galaxies with elevated star 
formation rates (SFRs) of a few - 100 Mq yr"^ (Berger et al. 2001, 2003a; Sokolov et al. 2001; 
Frail et al. 2002; Chary et al. 2002). To sustain such putative vigorous star formation, it is 
immediately clear that a large amount of molecular gas, a fuel of massive star formation, must 
exist in the host galaxies of GRBs. For example, a molecular gas mass of ~ IO^Mq is required 
in order to keep a SFR of 100 Mq y~^ for 10^ yr, which is a typical duration of starburst (e.g., 
Thornley et al. 2000, and references therein). It is also claimed that GRBs with undetected, 
or dark, optical afterglows ("dark bursts") probably occur in molecular clouds (Reichart, Price 
2002; but see Barnard et al. 2003); if it is the case, the surrounding interstellar medium (ISM) 
would be drastically affected by GRBs (Draine, Hao 2002). In fact, high column densities up 
to A^H ~ 10^2"^^ cm~^, consistent with the idea of GRBs within molecular clouds, are suggested 
based on an analysis of GRB afterglows at X-ray (Galama, Wijers 2001). 

Nevertheless, no observations of molecular gas have been reported toward GRB hosts, 
so far. This is basically due to the limited sensitivities of the current millimeter-wave facilities. 
Many of GRBs are at a cosmological distance {z> 1), but CO emission from galaxies beyond 
z = 1 is not very easy to observe, unless it is gravitationally amplified or extremely luminous 
(e.g., Barvainis et al. 2002). 

Here, we present a sensitive search for redshifted CO (J = 1 — 0) emission from the host 
galaxy of the luminous burst GRB 030329/SN 2003dh using the Nobeyama Millimeter Array 
(NMA). The continuum emission at millimeter wavelengths was also monitored simultaneously. 
GRB 030329 is one of the closest objects found among the luminous GRBs to date {z = 0.1685; 
Greiner et al. 2003; Caldwell et al. 2003), and therefore a suitable target for a deep search of 
molecular gas. 

Millimeter- wave continuum emission from afterglow is also of great interest. Because the 
spectral energy distribution (SED) of afterglow peaks at millimeter wavelengths on a timescale 
of a few days after a burst (e.g., Berger et al. 2000), flux monitoring at mm- wavelengths is very 
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important to determine the peak flux and frequency of the synchrotron spectrum as a function 
of time from the burst (e.g., Galama et al. 2000). Flux monitoring of afterglow at millimeter 
wavelengths has been made toward several bright GRBs, such as GRB 970508, GRB 980329, 
GRB 991208, GRB 000301C, GRB 010222, and GRB 030329 (e.g., Bremer et al. 1998; Taylor 
et al. 1998; Galama et al. 2000; Berger et al. 2000; Frail et al. 2002; Galama et al. 2003; Sheth 
et al. 2003; Kuno et al. 2004). 

Assuming a cosmology with Hq = 71 km s~^ Mpc~^, Qm = 0.27, and Q\ = 0.73, the 
luminosity distance of GRB 030329 is di^ = 805 Mpc, and the angular distance is = 589 
Mpc (1" corresponds to 2.86 kpc). 

2. Observations and Data Analysis 

We observed GRB 030329 using Nobeyama Millimeter Array (NMA) during the period 
from 2003 April 6 to 2003 May 30. NMA consists of six 10 m dishes equipped with cryogenically 
cooled receivers employing Superconductor-Insulator-Superconductor (SIS) mixers in double 
side band (DSB) operation. The most compact configuration (D array; baseline lengths ranged 
from 13 m to 82 m) was used for the observations, except for the first two days (C array 
for April 6 and 7; baseline lengths were from 26 m to 163 m). For 3 mm observations, the 
tracking frequency was set to 98.825 GHz to observe a redshifted CO(l-O) line (~ 98.65 GHz) 
at the upper side band (USB), whereas the lower side band (LSB), centered at 86.825 GHz, 
(and line-free channels of USB) was for the continuum. The center frequencies of USB and 
LSB for 2 mm observations were 146.969 GHz and 134.969 GHz, respectively. Each side band 
was separated by 90° phase switching of a reference signal. The Ultra Wide-Band Correlator 
(Okumura et al. 2000) was configured to cover a bandwidth of 1024 MHz per side-band with 
an 8 MHz resolution. 

A radio source, J1159-I-292, was observed every ~ 20 minutes for amplitude and phase 
calibrations, and the passband shape of the system was determined from observations of strong 
continuum sources, 3C 273 or J0423-013. We also observed J1058-I-015 several times during 
each observing run to check the consistency of the amplitude calibration. The flux density of 
Jl 159+292 was determined several times during the observing run; it ranged from 1.6 to 2.3 
Jy. In this paper, we adopted a constant flux value of 2.0 Jy throughout our observing run. 
This gave almost the same flux scale as the OVRO/BIMA results (Sheth et al. 2003; 1.95 Jy 
was adopted to the data). The overall error of the absolute flux scale was estimated to about 
± 20%. 

The raw visibilities were edited and calibrated using the NRO UVPROC-II package 
(Tsutsumi et al. 1997), and then USB and LSB data were concatenated to produce a final 
continuum image of each day using the NRAO AlPS task IMAGR. 

For a spectral line analysis, the continuum emission was subtracted first from the visi- 
bilities using the UVPROC-II task LCONT, if continuum was detected. Then, each data set 
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was combined into a single visibility set, and channel maps with a width of 32 MHz (97.2 km 
s~^) at an increment of 16 MHz (48.6 km s~^) were produced. 

3. Results 

3.1. Continuum Emission 

An unresolved 2 mm/3 mm continuum source was detected at the position of GRB 
030329 from 2003 April 6 to April 21. The observational log is provided in table 1. Some 
selected continuum images, showing a fading afterglow emission, are presented in figure 1, and 
the obtained millimeter- wave light curve is given in figure 2. 

The visibilities often suffered from large phase noises of atmospheric origin. To estimate 
the effect of the phase noise, we made images of the visibility calibrator source Jl 159+292 as 
well. We found that the decorrelation factors were typically ~ 0.95, which decreased to 0.90 
for the worst case. The obtained flux densities listed in table 1 are after a correction of the 
decorrelation. 

3.2. Search for Spectral Line Emission/ Absorption 

Neither significant emission nor an absorption feature was detected from individual spec- 
tra of each day. This is consistent with a report by Sheth et al. (2003). However, by summing 
up the data at 3 mm from 2003 April 10 to May 30, we found a possible emission feature 
from this integrated spectrum. Figure 3 displays the spectrum, and a velocity-integrated im- 
age of the possible emission feature. Note that the data on April 15 was not included in the 
final spectrum, because the weight of the data was very small due to its larger system noise 
temperature. 

The spectrum shows an emission feature with a peak flux density of 6.2 ± 3.8 mJy at 
98.62 GHz, corresponding to a redshift of zqq = 0.1689, if this is CO (J = 1 — 0) emission. The 
significance of the emission feature is low (~ 2a), but this feature appears on several successive 
channels, from zqo = 0.1682 to 0.1689. This is within the range of the reported GRB 030329 
redshifts determined from optical emission/absorption lines, 0.1685 (Greiner et al. 2003) to 
0.169 ± 0.001 (Price et al. 2003). 

The velocity width of this feature is 220 ± 49 km s~^. If we integrated the possible 
emission feature over the velocity width of 220 km s^^, we found a ~ 3(T source at the position 
of GRB 030329 (figure 3, right panel). The integrated CO flux density is 1.4 ±0.52 Jy km 
s~^. The position and redshift coincide well with those of GRB 030329, though we need further 
observations to confirm this possible emission-line feature. 
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4. Discussion 



4.I. Millimeter- Wave Light Curve: a Rapid Decay with Possible Plateaus? 

The overall appearance of our millimeter-wave light curve in table 1 and figure 2 is 
consistent with the OVRO/BIMA (Sheth et al. 2003) and NRO 45 m telescope (Kuno et al. 
2004) results; i.e., a rapid decay of millimeter- wave afterglow has been observed. Figure 4 shows 
a comparison of A = 3 mm measurements of GRB 030329 afterglow from NMA, OVRO/BIMA, 
and the 45 m telescope. A least-squares fit of all 3 mm observations during At = 9 to 23 gives 
a temporal decay index of asmm = —2.03 ±0.21, where millimeter fiux, S, decays as 5 oc t". The 
steep slope observed at millimeter-wavelengths, close to the decay indices after the jet break 
at At ~ 0.5 days at optical wavelengths (a = -1.97 ±0.12, Price et al. 2003; a = -1.81 ±0.04, 
Sato et al. 2003), suggests that the expansion during this period (At = 9 to 23) still remained 
relativistic (e.g.. Frail et al. 2004). 

In addition to the rapid decay, our data also implies that the decay was accompanied by 
possible "plateau" phases, or "bumps", with a time scale of ~ several days. One can see that 
the fiux densities in figure 2 are almost constant from At = 12 to 17 days (three days average 
is 17.6 mJy), whereas it drops suddenly at At = 18, and then the light curve becomes constant 
again till At = 23 (average is 9.0 mJy). The difference between these two plateau levels, a 
factor of 2, is significantly larger than the estimated overall fiux uncertainty, ~ 20%. In figure 
4, these possible plateaus/bumps are not very evident, however. This could be because the 
fiuxes in figure 4 were determined from different instruments, causing a significant error in the 
absolute cahbration. 

Note that the detection of fiuctuation at millimeter-wavelengths is important because 
the effect of interstellar scattering/ scintillation, which could be significant for radio light curves 
(e.g.. Frail et al. 2000), is mitigated at millimeter- wave light curves (Walker 1998). 

It is intriguing to see whether these possible bumps on the millimeter-wave light curve 
are achromatic or not, because simultaneous multi- wavelengths light curves provide a crucial 
clue to the origin of the bumps on the light curves (e.g., Nakar et al. 2003). Unfortunatelly, 
a detailed optical light curve compiled by Lipkin et al. (2004) lacks data points from At ~ 14 
to 21 days, where we find millimeter- wave bumps. Bloom et al.(2004) also presents late-phase 
(At > 10 days) optical/near-infrared light curves, yet their light curves also have gaps from At 
= 13 to 23 days. 

Some of the bumps or rebrightening events seen in optical/NIR bands are often ex- 
plained by a "refreshed shock" model (e.g., Granot et al. 2003), yet it remains unclear whether 
radio/millimeter- wave light curves can also be understood by this model (Berger et al. 2003b; 
Sheth et al. 2003; Frail et al. 2004). If the observed plateaus/bumps in the millimeter- wave 
light curve are caused by density fiuctuations of surrounding ISM (e.g., Wang, Loeb 2000), the 
spatial scale of the structure corresponding to the time scale of the observed bumps, several 
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days, is about ~ 10'^ AU, supposing a relativistic expansion. Small-scale structures of a dense 
molecular medium with such a spatial scale, a few 100 AU to 1000 AU, have been found in the 
Galactic star-forming regions (Kamazaki et al. 2001). 

4-2. CO Line from the Host Galaxy: a Very Gas/Dust Rich, Optically Faint Galaxy? 

Although the detection of the CO line from GRB 030329 should be confirmed by addi- 
tional/independent observations, if the emission-line feature in figure 3 is real, it gives valuable 
information about ISM in the host galaxy of GRB 030329. 

4-2.1. Molecular gas mass 

The integrated CO flux, 1.4 ±0.52 Jy km s~^, gives a coarse estimation of the molecular 
gas mass using a formula (Solomon et al. 1992), as follows: 



where aco is the CO luminosity to H2 molecular gas mass conversion factor, L'(^q is the CO 
line luminosity, S'co^'^ is the velocity-integrated CO flux, and di, is the luminosity distance. 
The CO line luminosity, L'qq, is 1.9 x 10^ [K km s^^ pc^]. The deduced mass reaches 5.5 x 10^ 
Mq, if a Galactic conversion factor, Xqo = 1-8 x 10^° [cm~^ (K km s"^)"^] (Dame et al. 2001) 
or aco = 2.89 in a unit of [Mq (K km s~^ pc^)~^], is applied. 

The molecular gas mass derived here is comparable to those in massive and luminous 
spiral galaxies (e.g., Nishiyama, Nakai 2001, and references therein), and significantly larger 
than the molecular gas masses in metal-poor dwarf galaxies. For instance, the total molecular 
gas mass of the Large Magellanic Cloud (LMC) is M(H2) = 5 x lO^M© (Fukui et al. 1999), 
adopting an Xco of 9 x 10^° [cm~^ (K km s~^)~^], determined from their own observations. 
Similarly, the molecular gas mass of the Small Magellanic Cloud (SMC) is M (H2) = 4 x IO^Mq 
by applying Xco = 2.5 x 10^^ [cm~^ (K km s~^)~^] (Mizuno et al. 2001), though this is a lower 
limit because their observation covered only a part of the SMC. 

Accordingly, the suggested molecular gas mass up to IO^Mq in GRB 030329 is sur- 
prisingly large if one recalls that the host of GRB 030329 is rather faint in optical bands; the 
absolute magnitude of the host galaxy is just comparable to that of SMC (Fruchter et al. 2003). 

One of the error sources of the deduced molecular gas mass is the CO luminosity to 
H2 molecular gas mass conversion factor. The molecular gas mass becomes even larger (by a 
factor of ~ 10) if we apply an Xco valid for low- metal galaxies (e.g., Arimoto et al. 1996). On 
the other hand, Xco factors derived in ultra-luminous infrared galaxies (ULlRGs) can reduce 
the gas mass by a factor of ~ 3 - 4 (Downes, Solomon 1998). In any case, supposing that the 



M(H2) = aco ■ ^co = "CO ■ 777 Sco^v ■ dl ■ i^J^ (1 + z) 




1.4 Jy km s 
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detection is real, a large amount of molecular gas, exceeding 10 Mq, must be present in the 
host galaxy of GRB 030329. 

4-2.2. Dynamical mass 

If the velocity width of the possible CO emission feature originated from rotation of the 
molecular gas, one can estimate an enclosed mass in the central region of the host galaxy by 
assuming the radius of a rotating gas disk and its inclination angle. The enclosed mass, often 
referred to as the dynamical mass Mdyn, is 



where r is the radius of the gas disk, V^tir) is an on-plane rotation velocity at a radius of r, 
AV is an observed line width, and i is the inclination of the disk {i = 0° is face-on). Here, we 
adopt the radius of the disk as r = 2" (or 5.7 kpc), an upper limit on the source size from our 
CO map, and an inclination angle of 45°, which is difficult to determine at this moment. In 
this case, the molecular gas mass shears about 17% of the total (dynamical) mass in the central 
kiloparsecs of the host galaxy. 

The major uncertainty of M^yn is inclination. If the disk is close to edge-on [i = 80° for 
instance), M^yn is reduced to be 1.6 x IO^^Mq. On the other hand, M^yn could exceed lO^^M© 
when the disk is closer to face-on {i < 23). Mgas/Mdyn is then ~ 1% to 10%, which is a typical 
range seen in nearby gas-rich spiral galaxies (e.g., Nishiyama, Nakai 2001). 

4.2.3. Molecular gas and star formation in the host of GRB 030329 

A gas-rich condition in the host galaxy of GRB 030329, suggested by our observations, 
implies the presence of an intense star-format ion/ starburst there because molecular gas masses 
roughly correlate with the rates of massive star formation (e.g.. Young et al. 1996). However, a 
SFR in the host of GRB 030329, determined from an extinction-corrected Ha luminosity, is 0.5 
Mq yrs~^ (Matheson et al. 2003), which is quite modest and not so enhanced compared with 
starburst galaxies (e.g., Kennicutt 1998). Here, we should note that the SFRs derived from the 
optical emission lines could underestimate the true rate of massive star formation if it is deeply 
embedded within dense molecular material. In fact, it is reported that the SFRs in some GRB 
hosts (such as GRB 980703, GRB 000210, GRB 000418, and GRB 010222), derived from optical 
observations, are much lower (by an order of magnitude) than those from radio/sub millimeter 
measurements, even after an extinction correction (Berger et al. 2001, 2003a; Frail et al. 2002; 
Gorosabel et al. 2003). This situation suggests that optical emission does not fully trace the 
whole massive star-forming regions in these GRB hosts; this could be the case in GRB 030329. 

For instance, it is likely that the unobscured SFR in the host galaxy of GRB 030329 




(1) 
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is 10-times larger than the SFR derived from optical data, as in the case of many GRB hosts 
(Berger et al. 2003). Then, the "true" SFR could be ~ 5 Mq yrs~^, and the expected radio (1.4 
GHz) and submillimeter (350 GHz) fluxes, originated from star formation in the host galaxy, 
would be ~ 0.3 mJy and ~ 0.1 mJy, respectively, according to fomulae by Carilli and Yun 
(1999). Such a flux at 1.4 GHz should be detectable using VLA when the afterglow emission 
has been disapperaed. 

Consequently, further multiwavelengths studies from radio to infrared, as well as addi- 
tional deeper CO observations, are indispensable to unveil the true properties of star formation 
and to conflrm the presence of a large amount of ISM, in the host galaxy of GRB 030329. 
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thank to NRO and all people who enabled us to carry out this target-of-opportunity program. 
We are grateful to the NRO/NMA staff for the operation and for continuous efforts to improve 
the NMA. We extend our gratitude to E. Nishihara for notification of GRB 030329, and to A. 
Endo for careful reading of the manuscript. KK was financially supported by The Mitsubushi 
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Fig. 1. Selected images of GRB 030329 afterglow at A = 2 mm and 3 mm to depict fading afterglow 
emission. The observed frequency, time after the burst, and the synthesized beam are indicated. The 
central cross in each panel marks the position of the GRB 030329 optical afterglow, aj2ooo = 10''44"50?03 
and (5j2ooo = +21'*31'18."15 (Uemura ct al. 2003). The contour levels arc —3 cr, 3 cr, 9 cr, and 15 a for the 
141 GHz image, and the contour intervals are 2 a for 93 GHz images. The tr of each day is listed in table 
1. 



100 



10 









































































1 












































* ^I" 

\ 






















































































1 1 

















10 30 50 70 

Time from the burst [days] 



Fig. 2. Time variation of A = 2 mm (square) and 3 mm (circle) continuum flux toward GRB 030329 
observed with the NMA. The upper limits arc 2 a. 
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Fig. 3. (Left) Observed spectrum at the 98 GHz band toward GRB 030329, produced by summing up 
the 3 mm data from 2003 April 10 to 2003 May 30. Tlie spectral resolution is 16 MHz or 48.6 km s"\ 
and the rms noise level is 3.8 mJy. A possible emission feature is seen at the redshift of GRB 030329 (z = 
0.1685). The velocity width of the possible emission feature is 220±49 km s^^. A horizontal bar indicates 
a velocity width of 200 km s^^ for reference. (Right) Total intensity map of a possible emission feature, 
integrated over a velocity width of 220 km s^^. The synthesized beam, shown in the right bottom corner, 
is 7."0 X 6"3 (position angle = 84°). The contour interval is 0.52 Jy beam^^ km s^^ (1 cr). The source 
coincides with the position of the GRB 030329 optical afterglow. 
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Fig. 4. Millimetcr-wavc (A 3 mm) light curve of the afterglow GRB 030329 taken with the NMA 
(circle; this work), OVRO/BIMA (diamond; Sheth et al. 2003) and NRO 45 m telescope (square; Kuno et 
al. 2004). A steep decline of the afterglow (5 oc i-2.03±o.2i^ indicated by a solid line), is evident. Insertion 
is a close-up view of the light curve for At = 9 to 30 days. 
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Table 1. Nobcyama Millimeter Array observations of GRB 030329. 
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(1) Observing date. (2) Beginning time of observation. (3) End time of observation. (4) Duration of the observation. Net on-sourcc 
integration time is typically about half of this. (5) Time since the burst (2003 March 29.48 UT). Middle of the observing run is 
taken here. (6) Observing frequency. An average of USB and LSB frequencies. "141 GHz" means an average of 146.969 ±0.5 GHz 
and 134.969 ±0.5 GHz bands, and "93 GHz" means an average of 98.824 ±0.5 GHz and 86.824 ±0.5 GHz bands, respectively (7) 
Flux density. (8) Noise level (1 a). (9) System noise temperature in DSB. (10) Atmospheric phase stability. "Bad" means a large 
fraction of the visibility data should be discarded due to large phase noise. 
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